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ABSTRACT 

High  resolution  electron  energs-  loss  spectroscopy  (HREELS)  measurements  base 
been  performed  on  Si(  1 1 1 )  surfaces  heavily  p-doped  by  the  decomposition  of  adsorbed 
decaborane  After  thermal  decomposition  of  the  decaborane  to  produce  B  atoms  on  the 
surface,  the  low  energy  electron  diffraction  pattern  shows  a  (VixVfIRJO'’  periodicity  due 
to  the  presence  of  1/3  monolayer  of  boron  in  the  second  complete  layer  The  HREELS 
data  exhibit  two  strong  features  ( 1 )  the  B-Si  dipole  mode  at  96  meV  and  <2)  a  broad 
surface  plasmon  mode  is  observed  at  - 1 00  meV  toss  energy  due  to  tbe  free  carriers  in  the 
region  below  the  B-reconslructed  surface  layer  We  have  investigated  the  electron  energy 
dependence  of  the  surface  plasmon  mode  in  order  to  determine  the  feasibility  of  using 
HREELS  to  determine  the  depth  profile  of  the  free  carriers  due  to  B  diffusion  into  the 
reition  100-1000  A  below  the  surface  Unexpectedly,  we  find  that  kinematic  factors  pla\ 
an  important  role  in  the  electron  energy  range  used.  1  5  -  28  eV.  and  thus  limit  the  degree 
of  quantitative  information  that  can  be  obtained  about  the  carrier  depth  profile  from 
HREELS  data 
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ABSTRACT 


High  resolution  electron  energy  loss  spectroscopy  (HREELS)  measurements  have 
been  performed  on  Si(l  1 1)  surfaces  heavily  p-doped  by  the  decomposition  of  adsorbed 
decaborane.  After  thermal  decomposition  of  the  decaborane  to  produce  B  atoms  on  the 
surface,  the  low  energy  electron  diffraction  pattern  shows  a  (V3xV3)R30°  periodicity  due 
to  the  presence  of  1/3  monolayer  of  boron  in  the  second  complete  layer.  The  HREELS 
data  exhibit  two  strong  features:  (1)  the  B-Si  dipole  mode  at  96  meV  and  (2)  a  broad 
surface  plasmon  mode  is  observed  at  -100  meV  loss  energy  due  to  the  free  carriers  in  the 
region  below  the  B-reconstructed  surface  layer.  We  have  investigated  the  electron  energy 
dependence  of  the  surface  plasmon  mode  in  order  to  determine  the  feasibility  of  using 
HREELS  to  determine  the  depth  profile  of  the  free  carriers  due  to  B  diffusion  into  the 
region  100-1000  A  below  the  surface.  Unexpectedly,  we  find  that  kinematic  factors  play 
an  important  role  in  the  electron  energy  range  used,  1.5  -  28  eV,  and  thus  limit  the  degree 
of  quantitative  information  that  can  be  obtained  about  the  carrier  depth  profile  from 
HREELS  data. 

1.  INTRODUCTION 

An  important  parameter  for  the  understanding  of  the  space-charge  region  near  the 
surface  of  a  semiconductor  is  the  doping  density  and  the  doping  concentration  profile  or 
dependence  on  depth.  Although  sputter-profile  analysis  such  as  Secondary  Ion  Mass 
Spectrometry  (SIMS)  can  give  the  atomic  concentration  profile,  direct  measurements  of 
the  carrier  concentration  profile  usually  require  contacts  and/or  interface  formation  with 
metals  or  other  semiconductors.  It  is  well  known  that  one  can  probe  the  carrier  density 
through  the  low  energy  free-carrier  plasmon  excitation  using  high  resolution  electron 
energy  loss  spetroscopy  (HREELS)  [1-5]  since  the  probing  depth  can  extend  to  -3000  A 
depending  upon  the  excitation  energy  and  the  energy  of  the  incident  electrons.  In  this 
paper  we  report  experiments  using  HREELS  to  probe  the  carrier-hole  plasmon  exciton  in 
the  space-charge  region  100-1000  A  below  the  surface  of  Si(l  1 1).  We  have  explored  the 
possibility  of  measuring  the  doping  profile  that  occurs  after  in  situ  boron  doping  of  the 
Si(l  1 1)  surface  by  thermal  diffusion  of  surface  boron  atoms  produced  by  thermal 
decomposition  of  adsorbed  decaborane  molecules  An  effective  carrier-hole  concentration 
on  the  order  of -10'^  cm'^  has  been  deduced  from  the  HREELS  measurements  for 
monolayer-level  boron  initially  deposited  on  the  surface.  The  surface  hole-plasmon  is 
highly  damped  due  to  the  low  mobility  of  carrier  holes  at  such  high  concentrations  and  due 
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to  the  strong  intravalence  band  transitions  resulting  from  the  small  spin-orbit  splitting  of 
the  valence  band  for  Si  (35  meV).  Under  these  conditions,  kinematic  effects  in  HREELS 
have  to  be  taken  into  account  in  order  to  give  a  quantitative  description  of  the  electron 
energy  loss  features.  Dipole  scattering  theory  in  conjunction  with  the  long-wavelength 
Thomas-Fermi  dielectric  function  gives  an  adequate  description  of  the  electron  energy-loss 
features.  However,  the  ability  for  HREELS  to  probe  different  regions  of  the  space-charge 
layer  below  the  surface  is  limited  by  kinematic  and  finite  spectrometer  acceptance  angle 
effects,  especially  at  low  incident  electron  energies. 

2.  EXPERIMENTAL 

The  HREELS  measurements  were  performed  in  a  UHV  system  with  a  typical  base 
pressure  of  IxlO'*®  mbar.  The  HREELS  measurements  were  carried  out  in  the  specular 
scattering  geometry  with  a  60°  incident  angle  at  electron  energies  1 .5  -  28  eV.  The 
HREELS  data  were  recorded  in  a  pulse-counting  mode  with  an  energy  increment  of  0.6 
meV  and  are  presented  without  any  post-acquisition  treatment.  The  boron  deposition  was 
carried  out  by  thermally  decomposing  decaborane  molecules  [6,7]  initially  adsorbed  on 
Si(l  1 1).  Subsequent  annealing  at  1200  -  1300  K  promoted  boron  diffusion  into  the  bulk 
of  Si(l  1 1).  The  details  of  this  study  are  reported  elsewhere  [7]. 

The  surface  structure  of  the  B-modified  (V3xV3)R30°  Si(l  1 1)  surface  is  shown  in 
Figure  1  a.  The  surface  Si  atoms  occupy  the  T4  adatom  sites  to  form  an  ordered 
(V3xV3)R30°  superstructure  on  top  of  the  (1x1)  bulk-terminated  Si  lattice.  The  B  atoms 
occupy  the  substitutional,  subsurface  S5  sites  in  the  first  Si  bilayer  directly  below  the 
surface  Si  adatoms  [8-10].  The  nominal  subsurface  B  concentration  is  0  33  per  (1x1)  unit 
cell  or  1/3  monolayer. 


The  existence  of  an  extensive  body  of  knowledge  on  boron  diffusion  in  Si  allows  a 
reliable  estimation  of  the  B  diffusion  length  in  Si(l  1 1)  under  our  experimental  conditions. 
For  the  B  concentration  of  interest  here,  the  diffusion  is  in  the  extrinsic  regime,  i  e.,  the 
diffusion  coefficient  is  concentration  dependent.  Based  on  the  annealing  temperature  and 
time,  an  empirical  diffusion  law  based  on  the  extrinsic  model  [11]  estimates  a  B  diffusion 
length  of  lO^-lCH  A.  This  distance  has  been  verified  by  secondary  ion  mass  spectroscopy 
(SIMS)  depth  profile  measurements  performed  separately  ex  situ  on  our  doped  Si(l  1 1) 
crystals  Since  the  typical  characteristic  wavevector  transfer  involved  in  dipole  scattering 
is  -10-2  A"'  in  the  present  HREELS  [1],  the  amplitude  of  the  excited  plasma  wave  decays  ’/ 
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exponentially  into  the  bulk  on  a  length  scale  significantly  smaller  than  the  B  diffusion 
depth.  Thus,  the  B-doped  region  below  the  Si(l  1 1)  surface  is  essentially  bulk-like  in 
term,?  of  its  dynamic  response  to  the  low  energy  electrons.  This  gives  rise  to  the  electron 
energy  loss  features  to  be  discussed  below. 


3.  RESULTS 

3.1  HREELS  Observation  of  the  B-Si  Vibrational  Mode  and  Surface  Hoie-plasmon 
Excitation 

After  boron  deposition  and  annealing  to  form  the  ordered  (V3xV3)R30°  surface,  a 
typical  HREEL  spectrum  taken  at  100  K  from  the  Si(l  1  l)-B(V3xV3)R30°  surface  is 
shown  in  Figure  lb.  There  are  two  very  distinct  energy  loss  features  in  the  spectrum.  The 
free-carrier  surface  plasmon  excitation  in  the  highly  B-doped  region  below  the  surface 
gives  rise  to  a  broad  energy-loss  feature  centered  at  (o^^.  Superimposed  on  this 
background  is  a  sharp  loss  peak  at  96  meV  which  belongs  to  a  B-Si  vibrational  mode 
(optical  phonon)  localized  in  the  surface  region.  In  an  earlier  work,  Rowe  and  coworkers 
[12]  have  characterized  this  mode  as  a  breathing  motion  of  B  in  the  S5  lattice  site  in  Si. 

Since  the  free-carrier  plasmon  frequency  is  proportional  to  (n  =  carrier 
density),  varying  the  doping  level  in  Si  is  expected  to  have  a  direct  effect  on  the  observed 
loss  peak  position,  tWpp,  due  to  the  surface  plasmon.  This  characteristic  is  verified 
qualitatively  in  the  HREEL  spectra  (Figure  2)  taken  at  a  constant  electron  energy,  Eq=6.  1 
eV,  on  selectively  B-doped  Si(l  1 1)  with  increasing  amounts  of  boron  as  judged  by  the 
Auger  intensity  ratios.  The  inset  of  Figure  2  shows  the  expected  dependence  on  boron 
density  (determined  by  Auger  intensity)  for  all  but  the  highest  exposure.  It  is  also  clear 
that  the  frequency  of  the  localized  B-Si  vibrational  mode  remains  unaffected  by  the 
variation  in  n.  Since  this  vibration  is  a  true  surface  mode  originating  from  the  vibration  of 
substitutional  B  atoms  against  the  surrounding  Si  neighbors,  the  vibrational  amplitudes  are 
localized  in  the  top  two  or  three  atomic  layers.  Therefore,  unlike  the  optical  phonons  (the 
Fuchs-Kleiwer  modes)  in  bulk  polar  semiconductors  such  as  GaAs  [13],  this  localized  B- 
Si  optical  phonon  is  not  expected  to  couple  to  the  carrier  surface  plasmon. 

3.2  Electron  Energy  Dependence  of  Hole-plasmon  Loss 
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Figure  3  shows  a  series  of  HREEL  spectra  taken  at  100  K  from  a  well-ordered 
(V3xV3)R30°  surface  at  various  primary  electron  energies,  Eq.  The  intensity  of  the 
localized  B-Si  vibration  decreases  monotonically  with  increasing  Eq.  This  behavior  is 
consistent  with  the  dipole  scattering  mechanism  which  predicts  an  Eq^  dependence  in  the 
loss  intensity  [1].  Another  prominant  characteristic  observed  from  these  HREEL  spectra 
is  that  the  loss  feature  due  to  surface  plasmon  excitation  experiences  a  pronounced  upshift 
in  peak  energy  with  increasing  Eq.  Since  the  peak-position  shift  exhibited  in  the  spectra 
appears  to  be  similar  to  what  was  observed  when  the  boron  concentration  was  changed,  it 
might  be  tempting  to  interpret  this  phenomenon  in  the  context  of  a  changing  carrier 
concentration,  n,  as  the  effective  probing  depth  of  HREELS  (to  be  defined  later)  varies 
with  £o  However,  it  has  been  recently  demonstrated  on  GaAs  [14,15]  that  the  kinematic 
factor  in  dipole  scattering  plays  an  important  role  for  surface  hole-plasmon  excitation. 
Thus  any  quantitative  interpretation  of  the  £o-dependent  spectra  obtained  here  has  to  take 
this  into  account.  In  order  to  fully  address  these  physical  issues  contained  in  our 
experimental  measurements,  we  next  present  a  model  dielectric  function  analysis  of  the 
electron  energy  loss  cross  section  based  on  the  well-accepted  dipole  scattering  mechanism. 
This  is  followed  by  a  comparison  of  numerically  simulated  energy  loss  spectra  with 
experimental  measurements. 

4.  Theoretical  Analysis 

4.1.  Model  Dielectric  Function 


We  adopt  a  long-wavelength,  local-dielectric-response  theory  upon  which  to  base 
our  model  calculation  and  consider  only  a  constant  carrier  concentration,  n,  for  simplicity. 
The  dielectric  response  of  the  free-carrier  holes  is  assumed  to  be  described  by  the 
frequency  and  wavevector  dependent  Thomas-Fermi  dielectric  function,  ej-p. 


Ajce  n 


(o. 


m[a}  +/a>r  -  D{k)\ 


co'  +  iQ)V  -  D{k) 


(1) 


with  D{k)  =  ,  where  Vp  is  the  Fermi  velocity,  oip  is  the  bulk  plasma  frequency  and 

is  the  high-frequency  dielectric  constant.  The  damping  factor,  T,  is  directly  related  to 
the  carrier-hole  mobility,  /r^,,  by 
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where  e  is  the  electronic  charge.  The  effective  mass  of  the  carrier-holes  is  taken  from  the 
contribution  of  both  light  and  heavy  holes  in  the  form  [14] 
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where  and  are  the  concentration  of  light  and  heavy  holes  in  Si. 


(3) 


In  the  limit  of  long- wavelength  and  zero-damping,  the  surface  plasmon  frequency, 
cOgp,  is  defined  as 


I  \l'2 

4;re*;j  _  <Up 
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(4) 


which  relates  co^  directly  to  . 

This  choice  of  dielectric  function  is  based  on  the  following  physical  considerations 
for  our  system:  (i)  the  A-range  involved  in  our  near-specular  HREELS  measurements  is  at 
least  one  order  of  magnitude  smaller  than  the  Fermi  wavevector,  Ap,  at  the  carrier  density 
(~10'^  cm'^)  under  study.  Thus  the  experiments  essentially  explore  the  dielectric  response 
of  the  system  in  the  long-wavelength  limit;  (ii)  at  such  a  high  doping  level,  the  Thomas- 
Fermi  screening  length  is  very  short  (a  few  A)  compared  to  the  depth  (~10^  A)  of  the  free- 
carrier  space-charge  region.  Thus  a  physical  condition  exists  to  apply  the  local  response 
formalism.  In  addition,  a  quantitative  comparison  between  the  local-response  and  a 
nonlocal  theory  [  1 6, 1 7]  under  accumulation  layer  conditions  indicates  a  very  promising 
agreement  between  the  two  at  long  wavelengths,  the  important  region  for  HREELS. 


4.2  Dipole  Scattering  Mechanism  and  Kinematic  Effects 


In  the  dipole  scattering  theory  which  describes  the  electron  energy  loss  process 
due  to  collective  excitations  at  a  solid  surface,  the  spectral  intensity  /(<y)  or  differential 
cross  section  (normalized  to  elastic  intensity)  is  represented  in  a  2-D  momentum  space 
integral  [1,3]  which  contains  a  product  term  of  the  kinematic  factor  and  the  loss-function 
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Im[-l/(e+l)]  integrated  over  the  parallel  momentum  transfer  range,  Qq,  defined  by  the 
spectrometer  acceptance  angle,  3c 


Ii(o)  = 


Qd'-Q 


q<Qc[w^^Q‘  +(<y-v,-0^] 


-[l  +  w(ty)]Im 


-1 


L£(<y,C)  +  lJ 


(5) 


where  ki  and  are  the  wavevectors  of  the  incident  and  scattered  electron  respectively  and 
«(<y)  is  the  Bose-Einstein  factor.  Inside  the  A-space  integral,  the  kinematic  factor  acts  as  a 
weighing  function  on  the  loss-function.  Im[-l/(e+l)],  at  each  parallel  momentum  transfer, 
Q.  If  we  are  in  the  parameter  regime  of  small  damping  where  the  loss-function  consists  of 
sharp  peaks  with  narrow  widths  (such  as  the  optical  phonons  in  an  ionic  crystal),  vary’ing 
the  kinematic  factor  only  modulates  the  intensities  of  the  energy  loss  peaks  in  the  final 
outcome  of  the  A-integral  in  /(<u)  (assuming  small  dispersion  of  e).  The  loss  peak 
positions  essentially  coincide  with  the  poles  in  Im[-l/(s+l)].  On  the  other  hand,  if  the 
loss-function  contains  a  broad  frequency  response  due  to  large  damping,  kinematic  effects 
will  then  become  important  in  determining  the  peak  position  in  the  differential  cross 
section,  /(<y).  This  type  of  loss-function  Im[-l/(8+l)]  is  encountered  in  the  current  study 
of  surface  hole-plasmon  excitation  in  the  heavily  B-doped  Si(l  1 1).  Again,  under  the 
conditions  of  negligible  dispersion  w^th  k,  the  energy  position  of  the  observed  loss  peak  in 
1(a))  differs  from  the  peak  position  of  Im[-l/(e+l)]  by  virtue  of  the  strong  (^-dependence 
of  the  kinematic  factor.  At  the  doping  levels  currently  achieved,  the  damping  width  of  the 
hole  plasmon,  T,  becomes  comparable  to  the  surface  plasmon  frequency,  o)^,  due  to 
strong  scattering  of  carrier-holes  by  the  ionized  dopants.  Under  this  circumstance,  the 
observed  loss  peak  position,  (Upp,  can  no  longer  be  simply  equated  to  o)gp  as  defined  by 
eqn.  (4),  and  a  full  evaluation  of  I(co)  is  generally  necessary  in  order  to  obtain  (y^p,  nf,  and 
r  from  the  measured  energy  loss  spectrum. 

Another  well-known  aspect  of  the  kinematic  effect  in  dipole  scattering  is  the 
selectivity  in  Q,  the  parallel  momentum  transfer,  with  varying  E^.  At  a  given  Eq  and  for  a 
given  energy  loss,  fico,  the  kinematic  factor  (or  the  dipole  lobe)  peaks  sharply  at  a 
particular  wave  vector  2~(^co/2£o)^o  Consequently,  a  region  of  the  space-charge  layer 
defined  by  an  effective  sampling  depth,  ',  is  selectively  probed  by  the  low  energy 

electron  [  1  ] .  It  has  been  suggested  [1,3,4,17]  that  this  property  can  be  utilized  as  a 
potentially  meaningful  experimental  tool  to  probe  the  dynamics  of  a  non-uniform  space- 
charge  layer  in  a  depth-resolved  fashion,  especially  in  the  near  surface  region  (1(P  -  10^ 

A).  However,  by  varying  the  incident  electron  energy,  Eq,  not  only  the  sampling  depth. 
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£4fr,  is  changed,  but  the  kinematic  factor  is  also  effectively  changed  as  well.  Under  the 
conditon  of  large  damping,  the  kinematic  effects  of  varying  Eq  should  also  be  considered 
in  interpreting  the  fg-dependent  plasmon  loss  spectra.  This  will  be  presented  below  with  a 
numerical  evaluation  of  eqn.(5)  using  parameters  best  representing  the  experimental 
conditions. 

5.  COMPARISON  BETWEEN  THEORY  AND  EXPERIMENT 

The  model  calculation  incorporating  both  the  dipole  scattering  kinematics  and  the 
model  dielectric  function  is  carried  out  numerically  by  assuming  a  circular  spectrometer 
aperture  and  a  fixed  value  of  3c=2.5°.  This  will  provide  the  basis  for  a  more  quantitative 
consideration  of  the  experimental  data  and  understanding  of  the  physics  contained  in  the 
HREELS  measurements. 

As  a  start,  we  apply  the  formalism  outlined  above  to  generate  a  theoretical  energy 
loss  spectrum.  The  carrier  density,  Nf,,  and  the  damping  factor,  T,  are  used  essentially  as 
fitting  parameters  to  the  experimental  spectrum.  Figure  4  shows  the  comparison  between 
the  measured  data  (points)  and  the  fit  (solid  line).  The  overall  quantitative  agreement 
between  theory  and  experiment  is  good.  It  should  be  pointed  out  that  the  theoretical 
curve  is  not  sensitive  to  9c  because  at  Eq=24  eV,  the  scattered  electron  angular 
distribution  is  very  narrow  and  the  dipole  lobe  peaks  at  9  <  0.2°  away  from  specular 
direction.  Virtually  all  the  scattered  electron  intensity  is  collected  by  the  spectrometer 
aperture,  9c  =2.5°,  used  in  the  calculation.  The  "best  fit'  yields  ti/j  =2.8x10*^  cm'^  and  r= 
74  meV,  which  corresponds  to  a  (ysp=102  meV.  We  treat  the  carrier  density  rtf, 
determined  in  this  fashion  as  an  "effective"  carrier  density  which  is  an  average  over  the 
distance  of  the  effective  electron  sampling  depth,  defined  earlier  At  Eg=14  eV,  this 
yields  c/efT*  ^  ^he  observed  plasmon  peak  energy.  The  damping  term,  T,  is  related 
to  the  carrier-hole  mobility,  /jf,,  through  eqn.  (2)  given  in  Section  4.1,  from  which  a  value 
/jf,-  34  cm^/V  sec  is  obtained  through  the  fitting.  This  value  agrees  well  with  the  reported 
hole  mobility  in  heavily  B-doped  Si  [18].  Also  plotted  in  Figure  4  (dashed  curve)  is  the  k- 
integrated  loss  function,  Im[-l/(e+l)],  which  is  used  to  illustrate  the  effect  in  the  absence 
of  the  kinematic  factor.  Its  behavior  shows  that  in  the  A-rangc  explored  in  HREELS, 
dispersion  in  has  only  a  negligibly  small  effect  on  the  energy-loss  peak  position.  The 
dominating  factor  for  shifting  the  peak  position  away  from  that  of  Im[-l/(e+l)]  is  the 
strong  G>-dependence  of  the  kinematic  factor.  It  is  clear  from  Figure  4  that  when  <ysp~r, 
the  full  spectrum  /(ry)  and  the  A-integrated  loss  function  peak  at  different  (U  values. 
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Incidentally,  neither  peak  position  corresponds  to  the  surface  plasmon  energy,  co^,  defined 
by  eqn.  (4).  For  precisely  this  reason,  the  determination  of  rtf,  from  the  observed  loss  peak 
position  is  no  longer  straightforward.  A  simplified  relation  linking  the  observed  loss  peak 
position,  tUpp,  and  the  surface  plasmon  frequency,  tu^p,  has  been  given  by  Meng  et  al.  [14] 
under  the  approximation  of  a  l/6>-like  behavior  of  the  kinematic  factor. 

Next,  we  use  the  same  fitting  parameters  to  generate  a  set  of  simulated  spectra 
with  varying  E^.  These  results  are  shown  in  Figure  5.  Varying  the  electron  energy,  Eq,  as 
noted  before,  effectively  varies  the  effective  electron  probing  depth,  and  the  kinematic 
factor.  It  can  be  seen  from  the  calculated  energy  loss  spectra  that  the  major  effect  is  a 
modulation  to  the  loss  intensities.  For  electron  energy  £o  ^  *oss  peak  position 

remains  essentially  unchanged.  Below  this  electron  energy,  a  widening  dipole  lobe  moving 
away  from  the  specular  direction  results  in  an  intensity  cut-off  on  the  high  energy-loss  side 
due  to  a  finite  spectrometer  acceptance  aperture,  This  high  energy-loss  cut-off  causes 
the  overall  weight  of  the  energy  loss  features  to  shift  toward  lower  energy-loss  side. 
Numerical  calculations  with  two  different  values  clearly  indicate  this  effect  as  shown  in 
Figure  5.  Even  with  a  constant  carrier  density  rtf,,  such  kinematic  effects  are  strong  under 
the  condition  of  large  damping  in  the  loss  function,  Im[-l/(8+l)]. 

Comparison  between  the  experimental  data  (Figure  3)  and  the  theoretical 
simulation  in  Figure  5  shows  reasonable  qualitative  agreement.  The  general  behavior  of 
the  loss  intensities  exhibited  in  the  experimental  spectra  are  fairly  well  reproduced  in  the 
simulation.  However,  the  experimental  spectra  appear  to  show  a  much  more  pronounced 
peak  shift  than  the  simple  uniform  model  dielectric  function  would  predict  with  a 
physically  reasonable  spectromet  angular  acceptance.  Such  a  downshift  in  plasmon  loss 
peak  energy  with  decreasing  Eq  can  be  clearly  seen  by  inspecting  the  measured  and 
simulated  spectra  in  the  electron  energy  range  6  -  28  eV  where  the  spectrometer 
acceptance  angle  effect  is  minimized.  This  suggests  that  there  is  also  a  significant 
contribution  from  a  non-uniform  carrier  concentration  profile  (such  as  a  depletion  layer)  in 
the  near  surface  region  (100  -  400  A)  that  enters  into  the  measured  energy  loss  spectra. 

But  the  full  analysis  is  difficult  due  to  the  small  spectrometer  acceptance  aperture  cut-off 
effect.  This  situation  exists  for  most  HREEL  spectrometers  in  operation  today. 

The  suggestion  of  a  non-uniform  carrier  concentration  over  the  HREELS  probing 
depth  (100  -  400  A)  is  further  substantiated  by  ejc  situ  dynamic  SIMS  analysis  performed 
on  this  Si(l  1 1)  crystal.  The  sputter-profile  measurement  verifies  the  existence  of  a  boron 
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concentration  gradient  which  is  characterized  by  a  sharp  rise  (about  one  order  of 
magnitude)  in  boron  concentration  over  the  initial  -300  A  in  depth  [19],  This  trend  is  in 
excellent  agreement  with  the  prediction  of  the  HREELS  analysis  presented  above. 

The  dependence  of  the  effective  electron  probing  depth,  on  the  incident 
electron  energy  has  been  discussed  in  the  past,  A  number  of  studies  (mostly  on  group  II - 
V  compound  semiconductors)  have  been  carried  out  [2-5,  13-15]  under  various  bulk 
doping  conditions.  Our  current  results  clearly  suggest  that  in  the  limit  of  small  damping 
(such  as  for  the  high-mobility  free  electrons  in  III-V  semiconductors),  it  is  feasible  to 
utilize  the  depth-probing  capability  of  HREELS  to  study  the  space-charge  region  r“ar  the 
surface  of  a  semiconductor  in  a  quantitative  and,  more  importantly,  non  'ntrusive  fash'on. 
However,  when  the  plasmon  damping,  T,  becomes  comparable  to  the  surface  plasmon 
frequency,  co^^,  kinematic  effects  in  HREELS,  as  discussed  above,  will  limit  the  extent  of 
quantitative  physical  information  that  can  be  extracted  from  the  observed  energy  loss 
spectra. 

Finally,  we  briefly  discuss  the  role  of  the  ordered  (V3xV3)R30°  subsurface  B  layer 
which  has  been  largely  ignored  in  previous  treatments.  As  a  p-type  dopant  in  bulk  Si,  B 
produce  charge  carriers  (holes)  in  the  Si  valence  band  via  thermal  exitation  of  valence 
electrons  into  the  B  acceptor  levels  near  the  edge  of  Si  valence  band.  However,  we  argue 
that  based  upon  the  unique  surface  electronic  structure  of  the  boron-modified  Si(l  1 1),  the 
ordered  1/3  ML  B  atoms  are  electrically  inactive  in  terms  of  fulfilling  their  usual  acceptor 
function.  Because  of  the  termination  of  the  bulk  crystal  structure  and  the  (V3xV3)R30° 
surface  reconstruction,  these  subsurface  B  atoms  no  longer  provide  similar  acceptor  levels 
near  the  edge  of  the  Si  valence  band  as  bulk  dopant  B  atoms  do  The  B-induced  band-gap 
surface  states  lie  1.3-1.8  eV  above  the  valence  band  maximum  (VBM)  [20]  and  remain 
unoccupied  at  100  K  while  the  electron  density  originally  associated  with  the  surface  Si 
adatoms  redistributes  toward  the  subsurface  B,  producing  an  occupied  surface  state  below 
VBM.  As  a  result,  no  charge  carrier  (hole)  is  produced  in  the  bulk  Si  valence  band  by 
these  specially-coordinated  subsurface  B  atoms.  In  the  frequency  range  of  interest  here, 
the  dielectric  response  of  this  near  surface  layer  only  comes  from  the  B-Si  optical  phonon 
plus  the  contributions  from  the  ion  cores. 
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6.  SUMMARY  AND  CONCLUSIONS 


We  have  used  HREELS  to  explore  the  dielectric  response  in  the  near  surface 
region  of  a  highly  B-doped  Si(l  1 1)  crystal.  The  main  conclusions  can  be  summarized  as 
follows: 

1 .  Following  the  thermal  diffusion  of  boron  in  the  Si(l  11),  surface  plasmon  excitation 
due  to  the  free-carrier  holes  has  been  observed  by  specular  HREELS  at  100  K  in  the 
electron  energy  range  1.5-28  eV.  The  hole-plasmon  is  highly  damped  due  to  the 
scattering  of  carriers  by  the  ionized  dopants  and  contributions  from  intravalence  band 
transitions. 

2.  The  dipole  scattering  theory  combined  with  a  Thomas-Fermi  model  dielectric  function 
has  been  shown  to  provide  a  proper  description  of  the  electron  energy  loss  cross 
section.  Kinematic  effects  have  to  be  taken  into  account  in  order  to  make  a 
quantitative  interpretation  of  experimental  electron  energy  loss  features  and  deduce  the 
surface  plasmon  frequency, 

3.  Evidence  for  a  non-uniform  carrier  concentration  in  the  near  surface  region  (<1000  A) 
has  been  suggested  by  comparison  of  HREELS  data  and  model  calculations. 

4.  At  very  low  electron  energies  (Eq  <  5  eV),  the  combination  of  strong  damping  and 
finite  spectrometer  acceptance  angle  cut-off  complicates  the  quantitative  interpretation 
of  the  experimental  HREELS  data. 
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Figure  Captions 


Figure  1 :  (a)  The  structure  of  the  Si(l  1  l)-B(V3xV3)R30°  surface  with  the 

(V3xV3)R30°  unit  cell  outlined.  The  B  atoms  occupy  the  subsurface,  substitutional  S5 
sites  in  the  first  Si  bilayer,  (b)  HREEL  spectrum  taken  from  a  well-ordered  Si(l  1 1)- 
B(V3xV3)R30°  surface  at  100  K  with  a  primary  electron  energy  Eq=3.2  eV.  The  peak 
position  corresponding  to  the  surface  hole-plasmon  loss,  rUpp,  is  at  63  meV.  The  loss  peak 
belonging  to  the  localized  B-Si  optical  phonon  mode  is  at  96  meV.  The  inset  shows 
schematically  the  composition  of  the  Si(l  1 1)  crystal  with  the  1/3  ML  B  layer  displayed  on 
an  expanded  scale. 


Figure  2:  HREELS  measurements  of  surface  hole-plasmon  excitation  on  selectively 

B-doped  Si(,l  1 1).  The  incident  electron  energy  is  6. 1  eV  at  a  60°  incident  angle.  The 
spectra  correspond  to  increasing  boron  concentration  in  the  near  surface  region  following 
sucessive  boron  deposition  and  annealing  cycles.  The  inset  shows  a  plot  of  (AE)^  versus 
the  B/Si  Auger  ratio. 


Figure  3:  HREELS  measurement  from  a  well-ordered  Si(l  1 1)-B  surface  as  a 

function  of  incident  energy,  Eq,  at  100  K.  The  dashed  line  marks  the  frequency  of  the  B-Si 
vibrational  mode,  v(B-Si). 


Figure  4:  Comparison  between  theory  and  experimental  HREELS  measurement. 

The  solid  theoretical  curve  is  based  on  the  Thomas-Fermi  model  dielectric  function,  Cyp, 
and  dipole  scattering  theory.  The  best  fit  to  the  e,:.perimental  spectrum  is  obtained  with 
«;,=2.8xl0’^  cm‘-^  and  a  damping  factor  r=  74  meV,  which  corresponds  a  surface  plasmon 
frequency  £t;5p=102  meV  and  a  loss  peak  position  (Upp=88  meV  The  dashed  curve  is  the 
Ac-integrated  loss  function  Im[-l/(s+l)].  A  Gaussian  function  with  1 1  meV  FWHM  has 
been  used  to  simulate  the  elastic  peak. 


Figure  5:  Calculated  energy  loss  spectra  I{co)  as  a  function  of  Eq  based  on  the  best-fit 

parameters  discussed  in  the  text.  The  solid  curves  are  for  a  spectrometer  acceptance  angle 
dc=2.5°,  the  dashed  ones  for  0° 
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